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bJO 1 Abstract 

We study limits of vacuum, isotropic universes in the full, effective, four-dimensional 

> : 

■<^j- . theory with higher derivatives. We show that all flat vacua as well as general curved 

rn ■ 

' ones are globally attracted by the standard, square root scaling solution at early times. 

Open vacua asymptote to horizon-free, Milne states in both directions while closed 

m . 

universes exhibit more complex logarithmic singularities, starting from initial data sets 



- 1—1 

X 



of a smaller dimension. 



^ ; 1 Introduction 

Vacuum states are not trivially obtainable for simple isotropic universes in general relativity, 
and one has to go beyond them to anisotropic, or more general inhomogeneous cosmologies 
for a vacuum to start making sense [I] . However, isotropic vacua are very common in effective 
theories with higher derivatives, see e.g., [2]-[Z]- Such classical vacua are usually thought of 
as acquiring a physical significance when viewed as possible low-energy manifestations of a 
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more fundamental superstring theory, although their treatment shows an intrinsic interest 
quite independently of the various quantum considerations. 

In this paper, we consider the asymptotic limits towards singularities of vacuum uni- 
verses coming from effective theories with higher derivatives. Such a study is related to the 
existence and stability of an inflationary stage at early times in these contexts, and also to 
the intriguing possibilities of solutions with no particle horizons. For flat vacua, we find the 
general asymptotic solution with an early-time singularity. This result is then extended to 
cover general curved vacuum isotropic solutions and we give the precise form of the attractor 
of all such universes with a past singularity. We also obtain special asymptotic states valid 
specifically for open or closed vacua starting from lower- dimensional initial data. 

The plan of this paper is as follows. In the next section, we derive the basic cosmological 
equations and show that they form an autonomous dynamical system in suitable variables. 
Then, in Sections 3-5, we apply asymptotic methods to study the early and late evolution of 
these isotropic cosmologies. In particular, we study separately the flat and curved subcases 
and show that there exist certain properties valid asymptotically irrespective of the influence 
of curvature, while other limits, coined here 'Milne states', have a very sensitive dependence 
on the sign of the constant curvature slices. 

2 The vacuum field 

We consider a vacuum, FRW universe with scale factor a(t) determined by the Robertson- 
Walker metric of the form 

g A = -dt 2 + a 2 g 3 . (2.1) 

Each slice is given the 3-metric 

93 = T ^ r2 dr 2 + r 2 g 2 , (2.2) 

k being the (constant) curvature normalized to take the three values 0,+l or —1 for the 
complete, simply connected, flat, closed or open space sections respectively, and the 2- 
dimensional sections are such that 

g 2 = d6 2 + sin 2 6d<p 2 . (2.3) 
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Our general higher order action is (we set 8irG = c = 1, and the sign conventions are those 

of m 

5 = 1 I C(R)dn g , C(R) = R + /3R 2 + -fRic 2 + SRiem 2 , (2.4) 
2 Jm 4 

where /3,j,6 are constants. We consider a family of metrics {g s : s G R}, and denote its 
compact variation by g pv = (dg/ds) s= Q. Since in four dimensions we have the Gauss-Bonnet 
identity, 

S GB = / (R 2 GB dn g y = 0, R 2 GB = R 2 - 4Ric 2 + Riem 2 , (2.5) 

J M 4 

in the derivation of the field equations through a g- variation of the action (12.41) . only terms 
up to Ric 2 will matter. Below we focus in the case where AA A is a homogeneous and isotropic 
universe with metric ( 12. ip . in which case we have a second useful identity, 

/ {{R 2 - 3Ric 2 )dfi g y = , (2.6) 
Jm 4 

which further enables us to include the contribution of the Ric 2 term into the coefficient of 



R 2 , altering only the arbitrary constants. Hence, the field equations read as follows: 

0, (2.7) 



2 y 6 



2RRT - ^R 2 g^ - 2(g' Mp g ua - g^g pa )V p V a R 



where we have set £ = 2(3/3 + 7 + 5). We note that because of the form of the coefficient 
£, some 'memory' of the original fully quadratic theory ( 12. 4p remains, and the final effective 
action leading to the field equations ( 12. 7ft is not equivalent to an R + (R 2 action with ( 
arbitrary (a use of the latter action would mean taking into account only the algebraic 
dependence of the quadratic curvature invariants). 

Eq. (12. 7p naturally splits into 00- and ^-components (i = 1,2,3), but only the 00- 
component of ( 12. 7p will be used below. Using the metric ( 12. ip . the field equation ( 12. 7p takes 
the form 



k + a 



;,2 



"d a da 2 d 2 d 4 , d 2 k 2 

2 — + 2 — - -- 3 - - 2fc - + - 

a 2 a A a 1 a 4 a 4 a 4 



0, (2. 



which after setting H = d/a for the Hubble expansion rate leads to our basic cosmological 
equation in the form 

1 H 2 k 1 k 2 1 1 k 1 / \ 

H=---MH + - 2 H--- i ---H-— 2 - (2.9) 



where now we have put e = £/6. Setting x = H, y = H and z = a 2 , Eq. (12. 9p can be 
written as an autonomous dynamical system of the general form 

x = f V Ac(x), x=(x,y,z), (2.10) 



that is 



x = y 

y = -3 xy + kxZ-— — — (2.11) 

2x 2x 12e 12ex 

z = —2xz 

equivalent to the vacuum, 3-dimensional vector field fyAC : K 3 — >■ M 3 : (x, y, z) >->■ f a;,rad(^, V, z) 
with 

y 2 k 2 z 2 x kz 



fvAc(x, y, z )=^y,--3xy + kxz- — - — -—, -2xz j . (2.12) 

This field completely describes the dynamical evolution of a vacuum, flat or curved, FRW 
universe in the gravity theory defined by the full quadratic action (12.41) . In the following, we 
shall assume that i^O, that is we consider only non-static universes. 

3 The unique flat vacuum 

When k — 0, the vacuum field ( 12.121) becomes 

(y 2 x \ 

y, 2x~ 3xy ~12e , ~ 2xZ )' ( ^ 3 ' 1 ' ) 

and the system (12. lip reads, 

x = y 

2 

y x 

V = 2x- 3xy -l2e (3 - 2) 
z = —2xz 

Our main interest below is to study the behaviour of the universe described by (13. ip . (13.21) 
near the initial singularity, taken here to lie at t — (it is really arbitrary, however, and we 
could have placed it at any to and used the variable r = t — to instead of t). Following the 
method of asymptotic splittings of Refs. [TU1 [HI H2], we find that of the 2 3 — 1 = 7 possible 
asymptotic decompositions that the field (13. ip possesses, there is only one that leads to a 
fully acceptable dominant balance, namely, 

f — f(°) _l_ f ( sub ) (O Q\ 

lO.VAC - r o,VAC + r 0,VAC \ 6 - 6 ) 
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with with dominant part 



f 0°VAc( X ) = ( V: ^ - 3x V' ~ 2xZ ) ' ( 3 ' 1 ) 



and sub dominant part 



f&A } c(x)= (0,-^,0). (3.5) 



We will eventually construct asymptotic series solutions which encode information about 
the leading order behaviour of all solutions as well as their generality (number of arbitrary 
constants) near the spacetime singularity at t = 0. For any given dominant asymptotic 
solution of the system (13. 2p . we call the pair (a, p) a dominant balance of the vector field 
fo.vAC) where a = (8,T],p) G C 3 are constants and p = (p,q,r) G Q 3 , and we look for a 
leading order behaviour of the form 

x(t) = ai p = (0t p ,rft q ,pt r ). (3.6) 

Such behaviours correspond to the possible asymptotic forms of the integral curves of the 
vacuum field fo,vAC, as we take it to a neighborhood of the singularity. Substituting the 
forms ( 13. 6ft into the dominant system (x,y,z)(t) = i*o°vac an< ^ s °l vm g the resulting nonlin- 
ear algebraic system to determine the dominant balance (a, p) as an exact, scale invariant 
solution, leads to the unique flat-vacuum balance £> 0) vac G C 3 x Q 3 , with 

i3 0iV Ac = (a,p)= (Q>"^p)> (-1,-2,-1)). (3.7) 

In particular, this means that the vector field f o° v ac ^ s a scale-invariant system, cf. [T0| [TTJ 
□2]. 

Further, we need to show that the term ( I3.5P in the basic decomposition (13. 3p of the 
flat- vacuum field ( 13. ip is itself weight-homogeneous with respect to the flat-vacuum balance 
(13. 7p for this splitting to be finally acceptable. For this we need to check that this candidate 
subdominant part is indeed subdominant. Using the balance £>o,vac defined by Eq. ( 13. 7p . 
we find that 

^^LfS-go (a) i> = (o,-£,o). (3.8) 

Taking the limit as t — > 0, this goes to zero asymptotically provided that the form fQ S y^ c (a) 
is different from zero. This happens only when e ^ 0, that is for all cases except when 
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3/3 + 7 + ^ = 0. We conclude that this basic decomposition is acceptable asymptotically in 
all higher order gravity theories when this constraint holds true. A physical example that 
is excluded in this analysis and consequently needs a separate treatment is the so-called 
conformally invariant Bach-Weyl gravity cf. [13]. Note that the same constraint appears in 
the stability analysis of the purely radiation universes in these theories, cf. Ref. [8]. 

We now move on to check our asymptotic solutions in terms of their consistency with 
the overall approximation scheme we are using (cf. [T0| [TT| IT2] for more details and proofs). 
We construct a series representation of the asymptotic solutions valid locally around the 
initial singularity, so that it is dominated by the dominant balance solutions we have built 
so far. The degree of generality of these final series solutions will depend on the number 
of arbitrary constants in them, and as explained in [10], the arbitrary constants of any 
(particular or general) solution first appear in those terms in the asymptotic series solution 
whose coefficients Cj have indices i = gs, where g is a non-negative /C-exponent and s denotes 
the least common multiple of the denominators of the set of all subdominant exponents and 
those of all the /C-exponents with positive real parts (in our case, s = 2). These exponents 
are numbers belonging to the spectrum of the so-called Kovalevskaya matrix given by 

/C = Df$ AC (a)-diag(p). (3.9) 

Hence, the /C-exponents depend on the dominant part of the vector field as well as the 
dominant balance. In the present case, the Kovalevskaya matrix is 

/ 1 1 0\ 



K. 



0,VAC 



(3.10) 



1 -1/2 
\ -2p / 

with spectrum 

spec(/C ,vAc) = {-1,0,3/2}. (3.11) 

The number of non-negative /C-exponents equals the number of arbitrary constants that 
appear in the series expansions. There is always the —1 exponent that corresponds to an 
arbitrary constant, the position of the singularity, and because the spec(/Co,vAc) i n our 
case possesses two non- negative eigenvalues, the balance £>o,vac indeed corresponds to the 
dominant behaviour of a general solution having the form of a formal series and valid locally 
around the initial singularity. To find it, we substitute the Fuchsian series expansions (no 
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constant first term, rational exponents) 

oo oo oo 

a:(t) = X)c li t5- 1 , y(t) = ^c 2 ^- 2 , z(t) =^c 3i ti- 1 , (3.12) 

i=0 i=0 i=0 

where Ci = 1/2, c 2 o = — 1/2, c 30 = p , in the original system (13. 2p and we are led to various 
recursion relations that determine the unknowns cu, c 2 j, c 3 j term by term. Further algebraic 
manipulations lead to the final series representation of the solution in the form: 

x{t) = \ t-' + cn t 1/2 --^ * + ••• (3.13) 
2 3oe 

The corresponding series expansion for y(t) is given by the first time derivative of the above 
expression, while the corresponding series expansion for z(t) is given by 

Z (t)=pt-^^t^ + ±-t + ... (3.14) 

Finally, we arrive at the following asymptotic form of the scale factor around the singularity: 

a(t) = a + 2 _^ t 2_^_ t 5/ 2 + 4«c|s f/2 + . . . (315) 
3 72c 9 

where a is a constant of integration and a~ 2 = p. 

As a final test for admission of this solution, we use Fredholm's alternative to be satisfied 
by any admissible solution. This leads to a compatibility condition for the positive eigenvalue 
3/2 and the associated eigenvector, 

v T ■ (k - 3 -I^j Cj = 0, (3.16) 

where I denotes the identity matrix, and we have to satisfy this at the j = 2 level. This 
gives the following orthogonality constraint, 



/ -|ci 3 + C 23 \ 



(2,1,-^ 



0. (3.17) 



ci 3 - 2c 23 

\ -2pCi 3 - |c 3 3 / 

It is immediate to check that this is indeed satisfied, thus leading to the conclusion that 
fl3.13p - (l3.14p corresponds to a valid asymptotic solution around the singularity. 

Our series solution f l3.13p - fl3.14p has three arbitrary constants, C13, p and another one 
corresponding to the arbitrary position of the singularity (taken here to be zero without loss 
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of generality), and is therefore a local expansion of the general solution around the initial 
singularity. Since the leading order coefficients are real, by a theorem of Goriely and Hyde, 
cf. [TT], we conclude that there is an open set of initial conditions for which the general 
solution blows up at the finite time (initial) singularity at t — 0. This proves the stability of 
our solution in the neighborhood of the singularity. 

We can see that the transformation C13 = 3c' 13 /2a and e = — k/6 in the series expansion 
(I3.15P leads to a form obtained by setting ( = in the series expansion found for the 
flat, radiation case of Ref. [Mj (cf. Eq. (21) in that reference, where the term 12£# 3 was 
mistakenly written as 24£# 3 there). 



4 Curved vacua 

As we have seen in the previous section when k = and we have a flat, vacuum FRW model 
in the fully quadratic theory of gravity defined by the action (12 .4p . the vector field fo,vAc 
has one admissible asymptotic solution near the initial singularity, namely, the form ( 13 . 1 3[) - 
(13. 141) . In this family, all flat vacua are asymptotically dominated (or 'attracted') at early 
times by the form a(t) ~ t 1//2 , thus proving the stability of this solution in the flat case. 

When k 7^ 0, and we have the situation of a vacuum, curved family of FRW universes, 
the vacuum field fyAc has more terms than those present in the flat case, namely, those 
that contain k in (12. 12ft . Below we shall use the suggestive notation f/^vAc instead of fyAc 
to signify that we are dealing with non-flat vacua. The field f^vAc (o r the basic system 
(12. lip ) can decompose precisely in 2 6 — 1 = 63 different ways. Of these 63 decompositions, 
there are only three that eventually lead to fully acceptable dominant balances, while the 
rest 60 decompositions fail to lead to an acceptable picture for various different asymptotic 
reasons. The only acceptable asymptotic splittings of the vector field ffc,vAc of the general 
form ffe 5 vAc — ffcvAC + ^feVAC' nave dominant parts 



-fe.VAC ' ^fe.VAC 

2x 



f fc°vACi = ( 2/> WZ ~ 3x V' ~ 2xz ) > ( 41 ) 



4°vac 2 = (^/' ~ 3x y + kxz > ~ 2xz ) , (4.2) 



f fe°VAc 3 = \ y '2^~ 3xy + kxz ~ ~2x~' ~ 2xz 1 ' ( 43 ) 
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while their subdominant parts are given respectively by the forms, 

f( s ^ r = ( 0, kxz - — - — , oV (4.4) 

^^(cg-^-^-^.o), (4,) 

«03- (0,-4-^,0). (4.6, 

The first decomposition (14. ip has identical dominant part as the flat splitting of the previous 
section, hence identical dominant balance, namely, its asymptotic balance is £>fc,vACi G C 3 x 
Q 3 , with 

B k ,y ACl = (a,p) = (Q,~,p) . (-1,-2,-1)) , (4-7) 

In particular, this means that the vector field f i°vACi * s a scale-invariant system. However, 
its subdominant part (14.41) is different, and we need to show that the higher order terms 
(I4.4p in the basic decomposition of the vacuum field are themselves weight-homogeneous 
with respect to the balance (14. 7ft for this to be an acceptable one. To prove this, we first 
split the subdominant part (14. 4p by writing 

= + f l 2 v ACl (x) + if VACl (x), ( 4 - 8 ) 

where 

Ca Ci (x) = (0,^,0), f£ VACl (x)= (°>-^-T^°)> f l!vA Cl (x)= ( '-^°)» 

(4.9) 

and using the balance £>fc,vACi defined by Eq. ( 14. 7p . we find that 
f(i) 



^fc,VACi( a ^ P / 
f(2) 



(4.10) 



f ^ VACl(atP) - fl 2 v ACl (a)t 2 = (0, (-kV - f, ) , (4.11) 

= fl 3 VACl (a)t 3 =(0,-^t 3 ,0). (4.12) 

Hence, taking the limit as f — > 0, we see that these forms go to zero asymptotically provided 
that ffcVACi( a )>* = 1,2,3 are all different from zero, which happens only when e 7^ 0, that 
is for all cases except when 3/3 + 7 + 5 = 0. Since the subdominant exponents 

q<P) = < g« = 1 < q® = 2 < q^ = 3, (4.13) 

9 



tP -i 



are ordered, we conclude that the subdominant part (I4.4p is weight-homogeneous as promised. 
Further, since the Kovalevskaya matrix and its spectrum are identical to the flat vacuum 
case, we arrive at the following asymptotic series representation for the decomposition ( 14. ip : 

<t) = lt- 1 -^ + c 13 t 1 /*-( , ^£ + ±-) « + ••-, (4.14) 



2 2 V 4 36e , 

while the corresponding series expansion for y(t) is given by the first time derivative of the 
above expression, and that for z(t) is given by 

z(i) . prl _V_^ ( v 2 + (Wl±M + _L) t + .... (4 . 15) 

For the scale factor, we find 

a(t) = at^-^t^ + ^t>-(^ + ^-) t^ + ..-, (4.16) 
2 3 \ 8 72e J 

where a is a constant of integration and a~ 2 = p. This series (I4.14p has three arbitrary 
constants, p, Ci 3 and a third one corresponding to the arbitrary position of the singularity, 
and is therefore a local expansion of the general solution around the initial singularity. The 
transformation C13 = 3c' 13 /2a and e = k/Q in the series expansion ( 14. 16ft leads to the form 
which is obtained by setting £ = in the series expansion found for the curved, radiation 
case, cf. Eq. (4.13) of [8]. In addition, by setting k = we are lead to the form (13. 13[) found 
for the flat vacuum. 

We note that because of the square root, limits can only be taken in the backward 
direction, t 4- 0, in the solution (I4.16p . another way of expressing the curious fact that this 
solution (along with Eq. (13. 15ft found in the previous Section) is only valid at early times 
and corresponds to a past singularity. 



5 Milne states 

We now move on to the analysis of the last two decompositions, namely, those with dominant 
parts ( 14. 2p and ( 14.3D . We show below that these lead to particular solutions for k = —1 and 
k = +1. In the case of open universes, k = —1, and the dominant parts take the forms 

f( °i,vAc 2 = (f > ~ 3x y ~ xz > ~ 2xz ) ' ( 5 - x ) 

f -ivAC 3 = (v, 7^ - 3xy ~ xz ~^ -2a*) , (5.2) 
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with subdominant parts given by 



{sub) I n V 



n(SU 



,VAC 2 



0, 



X z 

+ — ,0 , 



2x 2x 12e 12xe 



{{ sub) _( Q Q ) 
t-l,VACa-^, 12e + 12a . c »U;» 



respectively. These two forms lead to the same acceptable asymptotic balance 

S-i,VAC a , 3 = (a, P) = ((1, -1,1), (-1, -2, -2)) , 

while the structure of the /C-matrices is 

/ 1 1 \ 

£-l,VAC 2 — 



(5.3) 
(5.4) 

(5.5) 



2 -1 -1 
\ -2 / 



spec (AC. 



■1,VAC 2 



) = {-1,-1,2}, 



(5.6) 



and 



/C-i 



,VAC 3 



spec(/C_i,vAc 3 ) = {-1,-2,2}. 



(5.7) 



/ 1 1 \ 

2 -2 -2 

\ -2 y 

Since we are interested in the behaviour of solutions near singularities, we set the arbitrary 
constants corresponding to the negative eigenvalues equal to zero, and we are led to the 
following form for x(t), common for both decompositions, 



x(t) =t 1 + C12 t 



C12 - 18 e c\ 2 - : , 



60e 



(5.8) 



The corresponding series expansion for y(t) is given by the first time derivative of the above 
expression, while the corresponding series expansion for z(t) is given by 

'ci2 (42eci2 + 1 



Z (t) =t - C12 + 



120e 



t 2 + 



(5.9) 



Finally, we arrive at the following asymptotic form for the scale factor a(t) around the 
singularity: 

(5.10) 



q(i) = a t+^t 3 - a(Cl2 - l8e ^ t 5 + 
v ; 2 240e 



where ct = ±l as dictated by the definition z(t) = l/a(t) 2 . This solution has therefore 
two arbitrary constants, C\2 and a second one corresponding to the arbitrary position of the 
singularity (taken here to be zero without loss of generality), and is therefore a local expansion 
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of a particular solution around the singularity. Since the time singularity can be approached 
here from either the past or the future direction, we conclude that it represents a 2-parameter 
family of past, or future Milne states for these open vacua. This is also reminiscent of 
the Frenkel-Brecher horizonless solutions, cf. [15J, with the important difference that their 
solutions are matter-filled an possibly valid only in the past direction. 

On the other hand, when k — +1, the decomposition ( 14. 2 p does not lead to an acceptable 
dominant balance, but ( I4.3P does, namely, 



,(0) 



-1,VAC 3 



2x 



with subdominant part 



and we obtain 



-(sub) 
: +l,VAC 3 



3xy + xz , —2xz 

2x 



x z 
V2~e ~ 12xe'' 



£?+i,vac 3 = (a, p) = ((1, -1, 3) , (-1, -2, -2)) . 
The corresponding /C-matrix is 



(5.11) 

(5.12) 
(5.13) 



(£+i,vac 3 ) = {-1, -2V3, 2\/3}, 



(5.14) 



( 1 1 \ 

/C + i jV ac 3 = 10 -2 -2 , spec 

\ -6 ) 

and we expect particular solutions in this case with the given leading order, however, due to 
the irrational Kowalevskaya exponents the resulting series will contain logarithmic terms. 



6 Conclusion 

In this paper we have considered the possible singular behaviours and asymptotic limits of 
vacuum isotropic universes in the fully quadratic gravity theory which apart from the Einstein 
term contains terms proportional to a linear combination of R 2 , Ric 2 and Riem 2 . Taking into 
account various asymptotic conditions that have to hold in order to have admissible solutions, 
we are left with only three possible asymptotic decompositions of the vacuum vector field 
near the singular state. 

It turns out that a prominent role in the early asymptotic evolution of both flat and 
curved vacua in this theory is played by a scaling form that behaves as t 1 ! 2 near the initial 
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singularity. Using various asymptotic and geometric arguments, we were able to built a 
solution of the field equations in the form of a Fuchsian formal series expansion compatible 
with all other constraints, dominated asymptotically to leading order by this solution and 
having the correct number of arbitrary constants that makes it a general solution of the 
field equations. In this way, we conclude that this exact solution is an early time attractor 
of all homogeneous and isotropic vacua of the theory, thus proving stability against such 
'perturbations'. 

For open vacua, there is a 2-parameter family of Fuchsian solutions that is dominated 
asymptotically by the Milne form both for past and future singularities. In the case of 
closed models, we have logarithmic solutions coming from a manifold of initial conditions 
with smaller dimension than the full phase space but dominated asymptotically by the same 
a(t) ~ t form. 

It is instructive to also comment on our present results in connection with the recent 
results of [8] on the stability of radiation, curved universes for the same class of theories. 
We have shown that at early times both radiation and vacuum, flat or curved universes are 
past-attracted by the 'universal' t 1//2 asymptote and this is the most dominant feature in 
all these models. However, the behaviour of open vacua in these theories is more complex, 
for they allow universes that emerge from initial data sets of smaller dimension and valid 
for both early and late times. These universes asymptote to the Milne form during their 
early and late evolution toward singularities. Closed vacua, on the other hand, develop 
in time as more complex solutions that are characterized by logarithmic formal series, but 
asymptotically their leading order is described again by simple singularities similar to the 
open case treated here. 
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